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Foreword 


This  is  part  II  of  the  final  technical  report  for  Contract  No.  N62269- 
76-C-0378,  which  is  sponsored  by  the  Naval  Air  Development  Center,  Warminster, 
Pa.  The  work  was  performed  during  the  period  of  July  1,  1976  through 
December  30,  1977.  Mr.  Lee  W.  Cause  was  the  contract  monitor. 

The  contracted  study  is  under  the  title  "Certification  of  Composite 
Aircraft  Structures  under  Impact,  Fatigue  and  Environmental  Conditions"; 
parts  I  and  II  of  the  study  are  under  the  supervision  of  Dr.  P.C.  Chou, 
while  part  III  is  under  Dr.  A.S.D.  Wang,  both  of  Drexel  University. 

This  report  concerns  part  II  of  the  contract,  the  scale  effect  in 
fatigue  of  composite  materials.  It  is  a  self-contained  report  including 
definitions  of  all  nomenclature  used,  and  its  own  introduction  and  conclu¬ 
sions. 


The  authors  would  like  to  thank  Dr.  Edward  J.  McQuillen,  Dr.  James  L. 
Huang  and  Mr.  Lee  W.  Cause  for  the  frequent  technical  discussions.  The 
authors  would  also  like  to  thank  Mr.  James  Alper  and  Mr,  Dinh  Nguyen  who 
helped  conducting  the  experiments. 
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I.  Introduction 

It  is  known  that  the  endurance  limit  (fatigue  limit)  of  a  material  under 
fatigue  loading  decreases  as  the  size  of  the  specimen  increases  [1].  A  large 
amount  of  test  data  for  metals  are  available,  and  many  theories  have  been 
offered,  but  none  of  the  theories  are  satisfactory.  For  instance,  in  [2]  it 
was  found  that  for  an  Al-Cu  alloy  the  fatigue  limit  of  a  large  sheet  (230  mm 
width)  is  only  50%  of  that  of  a  small  sheet  (19  mm  width).  In  contrast,  the 
fatigue  limit  for  corresponding  mild  steel  specimens  showed  that  the  large 
sheet  has  85%  the  fatigue  limit  of  the  small  sheet.  Among  factors  mentioned 
that  may  contribute  to  the  scale  effect  include  maximum  stress  gradient,  heat¬ 
ing  of  large  specimen,  surface  stress,  etc.  No  statistical  theory  was  mentioned. 

Recently,  the  detrimental  effect  of  size  in  composite  material  has  been 
recognized  and  attempts  by  the  statistical  approach  have  been  proposed.  For 
instance,  in  Ref.  [3]  and  [4]  the  ratio  of  strength  of  two  specimens  was 
assumed  equal  to  the  ratio  of  their  volumes  under  stress.  This,  in  essence, 
is  considering  that  all  defects  or  micro-flaws  are  arranged  in  series,  and 
the  weakest  link  breaks  the  chain. 

The  concept  of  in-series  and  In-parallel  arrangement  in  statistical 
studies  is  not  new.  In  classical  reliability  analysis  the  series-parallel 
arrangements  are  often  used  for  engineering  systems,  see  for  instance  Shooman 
[5].  These  arrangements  have  also  been  assumed,  sometimes  implicitly,  in 
studying  the  strength  of  structures  and  materials,  see  [6].  It  seems  that 
they  have  not  been  used  to  study  the  scale  effect  of  fatigue  life. 

The  goal  of  the  present  study  is  to  obtain  basic  understanding  of  scale 
effect  in  the  fatigue  life  of  composite  material  structures.  This  basic 
understanding  may  be  used  to  develop  scaling  laws  such  that  the  service  life 
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of  a  full  size  structure  may  be  predicted  from  tests  of  small  size  models 
or  specimens . 

The  full  size  structure  is  considered  as  consisting  of  n  "elements". 

These  elements  may  be  statistically  in-series,  or  in-parallel,  or  a  coinbina- 
tion  of  these  two.  If  the  whole  structure  fails  when  any  one  of  the  elements 
falls,  then  these  elements  are  considered  to  be  in-series.  This  is  also 
known  as  the  weakest-link  theory.  On  the  other  hand,  if  after  one  element 
fails,  the  rest  of  the  elements  can  still  carry  the  total  load,  the  elements 
are  considered  in-parallel;  this  is  equivalent  to  a  bundle  of  loose  filaments, 
and  amenable  to  the  bundle  theory  analysis. 

If  the  strength,  or  fatigue  life,  of  the  element  has  a  small  statistical 
scatter,  or  large  Weibull  shape  parameter,  then  the  mean  strength  or  life  ^ 

of  the  n-element  structure  is  not  too  much  different  from  that  of  the  element. 

This  is  true  for  both  In-series  and  in-parallel  arrangements.  However,  if 
the  scatter  is  large,  which  is  the  case  for  fatigue  life  of  composites,  the 
n-element  structure  has  much  shorter  life  than  the  element. 

To  substantiate  this  in-series/in-parallel  model  in  relation  to  fatigue 
of  composites,  we  have  selected  as  the  basic  element  a  composite  specimen 
with  a  single  drilled  hole,  (Fig.  1).  By  arranging  a  number  of  holes  in  ^ 

different  arrays,  we  may  accomplish  in-series  and  in-parallel  models.  The 
present  paper  reports  the  results  of  the  basic  element  and  the  three  elements  , 

In-series  cases,  as  shown  in  Fig.  1.  This  drilled-hole  element  was  selected 
because  it  represents  the  rivet  or  bolt  holes  in  actual  structures,  and  it 
can  be  arranged  in  simple  in-series/in-parallel  models.  Fatigue  properties 
of  riveted  and  bolted  joints  in  composite  panels  are  reported  in  [7]  and  [8]. 
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Basic  Element  3  In-Series 

All  Dimension  in  Inch 
(1  Inch  =  25.4  mm) 

Properties 

Material:  Hercules  AS/3501-6  graphite  epoxy 
Layers :  12 

Layup:  1+  45/O2/+  ^5]^ 


Figure  1.  Specimen  Geometry  and  Properties 
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The  specimens  were  loaded  in  static  compression  and  constant  amplitude 
compression  fatigue.  The  compression  loading  is  used  because  it  is  known 
that  graphite/epoxy  composites  have  a  shorter  compression  fatigue  life  as 
compared  to  tension  fatigue,  see  Refs.  [9]  and  [10].  The  fatigue  tests  were 
at  a  constant  amplitude  of  =  0,  =  -2800  Ibf  (-12.50  kN)  (R  =  -«) . 

This  maximum  compression  stress  is  approximately  74%  of  the  mean  static  strength. 
Relevant  properties  and  dimensions  of  the  specimens  are  shown  in  Fig.  1. 

The  in-parallel  model  of  scale  effect  is  in  the  developement  stage. 

Pertinent  material  concerning  this  model  and  its  verification  is  contained 
in  Appendix  A. 
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II.  Equations  for  the  In-Serles  Model 

In  the  traditional  weakest  link  concept  in  material  strength,  the  basic 
"link"  is  a  flaw  in  the  material,  see  for  instance  Ref.  [11].  The  present 
in-serles  model  is  of  the  same  principle  as  the  weakest  link,  except  we  use 
a  structural  "element"  as  a  link.  This  element  can  be  of  a  size  in  the  same 
order  of  magnitude  of  the  total  structure.  It  is  well-known  that  the  weakest 
link  concept  can  be  described  in  mathematical  statistics  by  the  extreme  value 
theory,  see  for  instance.  Ref.  [12]  and  [13].  We  shall  summarize  the  results 
for  elements  of  Weibull  distribution. 

Consider  a  basic  element  of  Weibull  distribution,  with  a  density  function 

.a-1 


f(x)  = 


oCx-Xq) 


(1) 


and  a  cumulative  distribution  function  (CDF) 


Fj^(x)  =  P(X  £  x)  =  1  -  exp 


-'V  ) 


(2) 


where  x  =  strength  or  life,  Xq  =  position  parameter,  a  »  shape  parameter,  and 
3  =  scale  parameter.  Then  the  CDF  of  an  n-element  in-series  specimen  can  be 
shown  to  be 


F_  (x)  =  P(X  <  x)  =  1  -  [1  -  P(X  <  x)]“ 

The  density  function,  fjj(x) »  of  the  n-element  specimen  is  then 


uCx-Xq) 


a-1 


exp 


/  X-Xq  \ 

Ag/n^/V  J 


(3) 


(4) 


Comparing  Eq.  (4)  with  Eq.  (1),  we  conclude  that  the  shape  parameter  a,  and 

position  parameter  Xq,  of  the  n-element  in-series  model  are  identical  to  those 

of  the  basic  element;  the  scale  parameter  3  is  reduced  from  that  of  the  basic 

(1/h) 

element  by  a  factor  n 
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The  mean  strength  or  life  is  given  by 

"  er  (54-i)  +  (5) 

If  the  position  parameter  is  equal  to  zero  then  the  mean  of  the  n-element 
in-series,  y  is  related  to  the  basic  element  mean,  y  by 


The  decrease  of  the  mean  is  a  function  of  the  shape  parameter,  as  well  as 
the  number  of  elements,  n.  For  a  value  of  a  =  2,  the  strength  or  life  de¬ 
creases  a  great  deal;  for  n  =  100,  the  strength  of  the  in-series  combination 
is  only  10%  of  the  basic  element.  For  a  large  value  of  a,  however,  the 
decrease  in  strength  is  not  too  much.  For  instance,  for  n  =  100  and  a  =  30, 
y^  is  86%  of  y.  Eq.  (6)  is  plotted  in  Fig.  2. 

Since  for  composite  materials  the  shape  parameter  for  static  strength 
is  between  10  to  30,  and  for  fatigue  life  it  is  between  1  to  3,  it  can  be 
seen  that  the  size  effect  is  more  serious  in  fatigue  life  than  in  static 
strength,  if  the  elements  are  in-series.  In  other  words,  a  long  coupon  has 
almost  the  same  static  strength  as  a  short  element,  but  has  only  a  small 
fraction  of  the  fatigue  life  of  the  element. 
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Figure  2^  Ratio  of  Mean,  (of  a  Certain  Property)  of  n  Element  in  Series  oy^r  the  tj,ean  of  One  Element 
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III.  Static  and  Fatigue  Experiments 

To  verify  the  In-serles  model  a  number  of  static  and  fatigue  tests  were 
performed  using  both  the  basic  element  and  the  In-serles  specimen. 

A.  Specimens  All  specimens  used  for  this  study  were  cut  from  a  Hercules 
AS/3501-6  graphite /epoxy  composite  plate  prepared  at  the  Naval  Air  Development 
Center,  Warminster,  Pa.  The  plate  consisted  of  12  layers  following  a  layup 
of  [+  45/O2/+  ^5]g.  The  overall  measurements  of  the  plate  were  24  by  34 
Inches  (610  by  860  mm)  with  an  average  thickness  of  0.072  Inch  (1.83  mm). 

End  tabs  were  cut  from  plates  prepared  of  3M  Scotchply  1003  fiberglass /epoxy 

with  a  layup  of  [O2/9O/O2/9O/O2/9O/O2I •  The  thickness  of  the  plates  were 
0.080  Inch  (2.03  mm). 

The  basic  element  was  a  1  x  6  Inch  (25.4x152.4  mm)  strip  with  a  central 
1/4  Inch  (6.4  mm)  hole.  The  In-serles  specimen  had  the  same  overall  dimensions, 
but  with  three  1/4  inch  (6.4  mm)  holes.  End  tabs  of  length  1  1/2  inch 
(38.1  mm)  were  applied  to  the  ends  of  both  sides  of  each  specimen,  (See 
cutting  procedure  for  method  of  application).  The  dimensions  for  the' basic 
element  and  the  in-series  specimen  arie given  in  Fig.  1.  The  portion  between 
the  end  tabs  was  considered  as  the  test  section. 

The  plate  was  cut  Into  test  specimens  by  means  of  a  5  13/16  inch 
(147.6  mm)  diameter  diamond  saw  blade  mounted  on  a  Sundstrand  Rlgldmll. 

The  blade  was  run  at  715  rpm  with  a  table  feed  rate  of  1  inch  (25.4  mm)  per 
minute.  The  first  step  in  the  cutting  procedure  consisted  of  cutting  the 
plate  Into  6  1/2  x  10  Inch  (165.1x254.0  mm)  sections  such  that  the  0* 
fibers  were  parallel  to  the  shorter  edge.  Next  1  1/2  x  10  inch  (38.1x254.0  mm) 
end  tab  strips  of  0.080  Inch  (2.03  mm)  thickness  were  bonded  to  both  sides 
of  the  sections,  1/4  Inch  (6.4  mm)  from  each  of  the  10  inch  (254.0  mm)  edges. 
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The  bonding  agent  used  was  Ecobond  51  conbined  with  Catalyst  #9,  both 
manufactured  by  Emerson  and  Cuming  Inc.  These  sections  were  then  squared 
off  and  cut  down  to  6  x  10  inch  (152.4x254.0  mm)  by  taking  off  1/4  inch 
(6.4  mm)  along  each  of  the  10  inch  (254.0  mm)  edges.  Next  the  sections 
were  cut  into  6x1  inch  (152.4x25.4  mm)  strips.  Each  section  yielded  9 
of  these  strips. 

Specimen  holes  were  drilled  with  straight  fluted  solid  carbide  drill 
bits  manufactured  by  Cleveland  Twist  Drill  Co.  First  a  pilot  hole  was 
made  in  the  specimen  with  a  1/8  inch  (3.2  mm)  bit.  This  was  enlarged  with  a 
1/4  inch  (6.4  mm)  bit.  The  bits  show  no  noticeable  wear  after  drilling 
in  excess  of  200  holes. 

After  cutting  and  drilling  the  completed  specimens  were  stored  at 
room  temperature  and  humidity  until  testing. 

B.  Testing  Apparatus  All  testing  was  performed  on  an  Instron  Model  1230 
Dynamic  Test  Machine  equipped  with  a  Model  836  remote  control  panel,  a  Model 
602  load  and  stroke  controller  panel,  and  a  Model  860A  function  generator. 
This  testing  machine  has  the  capacity  of  +  10,000  Ibf  (+  40  kN)  load  which 
can  be  applied  in  the  form  of  a  sine  input,  triangular  input,  square  input, 
or  ramp  input.  The  maximum  reliable  frequency  for  operating  this  machine 
is  30  Hz.  The  Instron  Wedge  Action  Gripping  Jaws  were  used  in  conjunction 
with  our  own  gripping  and  anti-buckling  guide  system  (to  be  described  later) . 
For  static  tests  a.  Hewlett  Packard  Model  7045A  X-Y  Recorder  was  used. 
Additional  equipment  included  a  Tektronix  7623A  Storage  Oscilloscope  and  a 
Keithley  Instruments  160  Digital  Multimeter. 
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A  method  for  gripping  the  specimen  and  for  preventing  buckling  during 
loading  was  designed,  and  a  sectioned  assembly  of  this  apparatus  is  shown  in 
FiS*  As  can  be  seen  the  bottom  and  top  grips  are  held  in  place  by  the 
Instron  Jaws.  The  gaps  between  the  jaws  and  the  grips  are  filled  with  shim 
stock.  The  top  grip  is  stationary  and  the  two  vertical  runners  are  bolted 
lb.  The  bottom  grip  is  free  to  ride  up  and  down  in  these  runners.  The 
applied  load  is  transmitted  by  the  actxiator  arm  of  the  Instron  through  the 
bottom  Instron  Jaw  to  the  bottom  grip.  The  frictional  force  between  the 
bottom  grip  and  the  vertical  runners  was  found  to  be  negligible  (less  than 
10  Ibf  or  *t4  N) .  The  anti— buckling  guide  also  fits  in  the  runners  and  is  free 
to  move  up  and  down.  The  runners  Insure  correct  alignment  of  the  top  and 
bottom  grips  and  the  anti-buckling  guide.  Hence  the  test  specimen  will  have 
no  initial  crookedness.  The  clearance  between  the  anti-buckling  guide  faces 
and  the  specimen  surfaces  of  course  varies  from  specimen  to  specimen  but 
on  the  average  is  0.003  inch,  (0.08  mm).  The  section  view  shows  that  the 
specimen’s  ends  are  in  contact  with  the  grip  floors.  This  insures  that  the 
load  is  applied  on  the  ends  of  and  along  the  axis  of  the  test  specimen. 

The  figure  also  indicates  that  the  rear  gripping  plates  are  held  in  place 
by  1/4  20  Allen  screws  Inserted  through  the  front  grips.  The  rear  anti- 
buckling  guide  plate  is  held  to  the  front  anti— buckling  guide  by  four  screws . 

Detail  drawings  of  the  grips  and  anti-buckling  guide  are  presented  in 
Appendix  B. 

Tsstlng  Procedures  All  tests  were  run  at  room  temperature  with  ambient 
humidity.  In  order  to  reduce  friction  between  the  specimen  and  the  antl- 
buckllng  guide,  a  thin  coating  of  Molykote  44  grease,  a  heat  stable  silicone 
lubricant,  was  applied  to  the  contact  surfaces  of  the  guides. 
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Figure  3.  Sectioned  Assembly,  Specimen  Gripping  Anti-Buckllng  Device 
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For  static  tests  the  specimen  was  secured  In  the  grips  and  the  back 
of  the  anti-buckling  guide  was  screwed  in  position.  A  compressive  load 
of  10  Ibf  (40  N)  was  applied  to  the  specimen  and  then  the  screws  holding  the 
backs  of  the  grips  and  antl-buckllng  guide  were  torqued  tight.  The  load 
was  then  released.  Next,  operating  in  the  stroke  mode  and  using  the  ramp  ■ 
input  a  compressive  load  was  applied  to  failure.  The  failure  load  was  re^ 
corded  on  the  X-Y  plotter. 

For  fatigue  tests,  the  specimen  was  secured  in  the  grips  and  the  back 
of  the  anti-buckling  guide  was  screwed  in  place  (hand  tight) .  Small 
pieces  of  styrofoam  were  placed  between  the  guide  and  the  grips  in  order  to 
damp  out  any  vibration.  Operating  in  the  load  mode  a  compressive  load  was 
gradually  applied  up  to  100  Ibf  (440  N) .  Then  the  screws  holding  the  grips 
and  guide  in  place  were  torqued  tight.  The  load  was  then  Increased  slowly 
until  the  desired  mean  load  was  obtained.  In  the  present  test,  the  mean  load 
is -1400  Ibf  (-6.23  kN).  A  sine  wave  cycling  of  the  applied  load  with  an 
amplitude  of  1400  Ibf  (6.23  kN)  was  then  started.  For  the  first  10  cycles  a 
frequency  of  0.01  Hz  was  used.  During  this  time  the  loading  was  rechecked 
with  the  use  of  the  digital  multimeter  and  oscilloscope  to  make  sure  that 
the  specimen  did  not  go  into  tension  and  the  load  was  cycling  between  0  and 
—2800  Ibf  (-12.50  kN) .  The  frequency  was  then  gradually  and  slowly  increased 
to  5  Hz.  Petlodlcally  the  loading  was  checked  with  the  oscilloscope  and 
screws  retorqued  as  necessary.  The  specimen  was  run  to  failure  and  fatigue 
life  was  recorded  in  cycles  by  the  control  panel  counter. 
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IV.  Experimental  Results 

This  section  presents  the  static  and  fatigue  test  results  for  the 
basic  element,  the  three-ln-series  specimen,  and  several  other  arrangements. 
All  specimens  failed  through  a  drilled  hole.  For  each  test,  the  ultimate 
load,  or  the  fatigue  life  of  each  specimen  is  given  in  tabular  form.  The 
sample  means  and  sample  standard  deviations  are  also  calculated  from  the 


following  formulas 


Sample  mean 


1  " 

=  ~  I 

”  i=l  ^ 


Sample  standard  deviation  *  s 


where  n  is  the  total  number  of  specimens  tested  in  each  case. 


Table  I 

Static  Strength  and  Fatigue  Life  of  the  Basic  Element, 
Graphite/Epojsy  [+  45/O2/+  45]g 

Ultimate  Static  Compression  Load,  Ibf  (kN),  (12  specimens) 


2950  (13.12)  3600  (16.01)  3900  (17.35)  4350  (19.35) 

I  3100  (13.79)  3630  (16.15)  3900  (17.35)  4600  (20.46) 

3250  (14.46)  3770  (16.77)  3940  (17.55)  4650  (20.68) 

X  =  3800  (16.90) 
s  =  550  (  2.45) 

b.  Compression  Fatigue  Life,  Cycles,  Max.  Compression  Load  =  74%  Basic 
Element  Mean  Static  Strength,  R  =  Frequency  -  5  Hz  (10  specimens) 

17230  26230  78750  107880  127540 

20010  49460  83300  112120  128780 

X  =  75130 
s  =  44200 
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Table  la  presents  the  static  strength  measured  in  Ibf  (kN)  for  the  basic 
elements.  The  sample  mean  and  sample  standard  deviation  are  3800  and  550  Ibf 
(16.90  and  2.45  kN)  respectively.  Table  Ib  gives  the  fatigue  life  measured 
in  cycles  for  the  basic  elements.  The  maximum  compressive  load  to  which  each 
element  was  subjected  was  2800  Ibf  (74%  of  mean  static)  while  the  minimum, 
load  was  0.  The  mean  and  standard  deviation  of  the  basic  element  fatigue  life 
were  75,130  and  44,200  cycles  respectively. 


Table  II 


Static  Strength  and  Fatigue  Life  of  the  Three-In-Serles 

Specimen,  Graphlte/Epoxy  [+  45/0_/+  45] 

*•  s 


a.  Ultimate  Static  Compression  Load,  Ibf  (kN),  (10  specimens) 

3000  (13.34)T*  3260  (14.50)T  3450  (15.35)B  3660  (16.28)B 

3150  (14.01)E  3400  (15.12)B  3600  (16.01)B 

3150  (14.01)T  3400  (15.12)M  3600  (16.01)T 

X  =  3370  (14.99) 

s  =,  220  (  1.00) 

b.  Compression  Fatigue  Life,  Cycles,  Max.  Compression  Load  =  74%  Basic 
Element  Mean  Static  Strength,  R  “  -«>,  Frequency  =  5  Hz  (9  specimens) 

2090  M  15270  B  19990  M  29950  M  47690  B 

5040  T  16090  B  24700  B  32740  B 

X  =  21510 

s  =  14200 

*  Failed  hole  location,  T  =  top,  M  =  middle,  B  =  bottom. 


Table  Ila  shows  the  static  strength  measured  in  Ibf  (kN)  for  the  three-in- 
serles  specimens.  Also  the  position  of  the  failure  is  given.  The  sample 
mean  is  3370  Ibf  (14.99  kN)  and  the  standard  deviation  is  220  Ibf  (1.00  kN). 
Table  lib  presents  the  fatigue  life  in  cycles  and  the  position  of  the  failed 
hole  for  each  of  the  three-in-series  specimens  subjected  to  the  same  loading 
conditions  as  the  basic  element  fatigue  specimens.  Here,  x  =  21,510  cycles, 
and  8  “  14,200  cycles. 
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We  see  that  in  the  static  case  the  mean  strength  of  three-*in-series 
specimen  is  89%  of  the  basic  element  mean  strength.  However,  the  mean  life 
of  the  three-in-series  specimen  is  29%  of  the  basic  element  mean  life.  It 
is  evident  that  the  in-series  configuration  greatly  reduces  the  fatigue  life, 
but  has  a  much  less  effect  on  the  static  strength. 

Some  additional  tests  were  performed  on  specimens  with  different 
geometries  than  described  above.  This  data  is  compiled  in  Appendix  C. 
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V.  Parameter  Estimations 
A.  Interval  Estimation  for  the  Mean 

According  to  the  in-series  model,  the  life  and  strength  of  the  three- 
in-series  specimen  should  be  lower  than  those  of  the  basic  element.  From 
the  data  in  the  previous  section,  it  can  be  seen  that  indeed  the  sample 
mean  life  and  mean  strength  for  the  three-in-series  are  both  lower  than  the 
corresponding  basic  element  values.  The  sample  mean  values,  however,  are 
random  variables,  dependent  on  the  sample  size,  and  by  themselves  have  less 
statistical  significance.  In  order  to  establish  a  confidence  level  on  the 
hypothesis  that  the  three-in-series  and  basic  element  specimens  were  from 
two  different  populations  and  that  the  three-ln-series  case  has  a  lower  mean 
value,  we  shall  make  interval  estimations  for  the  means. 

We  shall  assume  that  the  random  variable 


has  the  Student-t  distribution  with  n-1  degree  of  freedom.  Strictly  speak¬ 
ing,  this  is  true  only  when  the  population  is  normal.  For  practical  pur¬ 
poses,  it  is  a  good  enough  approximation  for  non-normal  distributions,  [14]. 
Based  on  this  assumption,  we  have  estimated  the  confidence  Intervals  for  all 
four  cases,  and  listed  the  results  in  Table  III. 


■m 
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Table  III 


Confidence  Intervals  for  Basic  and  Three- In-Series  Specimens 
Fatigue  Life 


Specimen 

Confidence  Interval 
for  Mean  Life,  Cycles 

%  Confidence 

Basic 

[43510, 

106750] 

95% 

3- in-series 

[10600, 

32420] 

95% 

Static  Strength 

Specimen 

Confidence 

Intervals 

%  Confidence 

for  Mean  Strength 

Ibf 

kN 

Basic 

[3520,  4090] 

[15.7,  18.2] 

90% 

[3460,  4150] 

[15.4,  18.5] 

95% 

Three- in-s  er ies 

[3240,  3500] 

[14.4,  15.6] 

90% 

[3200,  3530]  [14.3,  15.7]  95% 


As  can  be  seen  the  95%  confidence  Intervals  of  the  population  mean  lives  of 
the  basic  element  and  three-ln-serles  specimen  are  disjoint,  and  the  forn^r 
larger  than  the  latter.  This  implies  that  these  sample  means  come  from 
separate  populations  and  the  basic  element  has  longer  life.  Table  III  also 
lists  the  90%  and  95%  confidence  Intervals  for  the  basic  and  three-in- 
series  mean  strengths.  Here  it  is  seen  that  the  95%  confidence  intervals 
overlap  whereas  the  90%  intervals  do  not.  It  can  be  argued,  although  not 
as  strongly  as  in  the  fatigue  case,  that  the  strength  results  of  the  basic 
and  three- in-serles  elements  do  indeed  come  from  different  populations. 

B.  Point  Estimation,  Weibull  Distribution 

It  is  well-known  that  the  strength  and  fatigue  life  of  materials  can 
be  best  characterized  by  the  log-normal  or  the  Weibull  distributions.  In 
this  study,  the  two  parameter  Weibull  distribution  is  used  to  represent 

our  results.  The  cumulative  distribution  ftmction  (CDF)  has  the  form, 
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F(x)  =  P(X  <  x)  =  1  -  exph  ] 


(10) 


Having  collected  sample  points  x,  the  parameters  can  be  estimated 
in  many  ways,  for  Instance,  method  of  moments,  maximum  likelihood,  and 
regression.  In  this  study  the  method  of  linear  regression  was  used  to  esti¬ 
mate  a  and  3.  This  technique  requires  an  assignment  of  numerical  values 
of  (called  the  rank)  for  each  x^^.  There  are  many  ways  of  making  this 
assignment  but  in  this  study  the  median  rank  method  was  used.  Here  the 
sample  points  are  ordered  from  lowest  to  highest  and  then  P^  is  assigned  by 
the  approximate  median  rank  formula 


p  -  1  -  0-3  .  _  ^ 

i  n  +  0.4*  “ 


(11) 


where  n  is  the  total  number  of  data  points  comprising  the  sample. 

In  using  linear  regression,  let  us  transform  Eq.  (10),  by  taking  the 
logarithm  twice,  into  the  form 


y'  =  ax'  -  a  Jin  3 


(12) 


where 

X*  =  Jin  X  (13) 

•  <“> 

At  X  -  x^,  the  corresponding  values  of  x^  and  y^  can  be  evaluated  from 
Eqs.  (13)  and  (14).  The  quantity  y^  is  associated  with  the  theoretical  dis¬ 
tribution  (10),  evaluated  at  x^.  The  corresponding  value  from  the  approximate 
rank  is  given  by 
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The  difference  which  represents  the  error,  is  squared,  and 


summed  over  i  or, 


I 

i=l 

n 

“I  (y4  -  ax'  -  a  In  3)‘ 

i=l 


This  sum  of  the  square  of  errors  is  then  minimized,  by  forming  the  two 
2  2 

equations  3  A  /3a  =  0  and  3A  /33  =  0.  Solving  these  two  equations  yields 
the  two  unknown  parameters  as 


n  n 


^  ^i  ^i  ■  n  ^  *1  ^i 


i«l  i«l 
n  v2 


3  =  exp 


_  n  n 

a  I  -  I  y^ 

i=l  i=l 


Table  IV  gives  the  values  of  a  and  3  for  the  four  cases.  Values  of 
the  mean  y  end  coefficient  of  variation  are  also  listed  and  were  obtained 


from  the  formulas 


y  =  3r 


c.o.v.  = 
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Tabic  IV 

Estimated  Weibull  Parameters  for 
Basic  Element  and  Three-In-Series  Specimens 

a.  Fatigue-Life,  cycles 


Specimen 

a 

a 

cycles 

y 

cycles 

C.O.V. 

Basic 

Element 

1.36 

87840 

80400 

.74 

Three- In-serles 

1.09 

25390 

24700 

.93 

b.  Static  Strength 

Specimen 

,  Ibf  (kN) 
a 

3 

Ibf  (kN) 

y 

Ibf  (kN) 

C.O.V. 

Basic 

Element 

7.6 

4040 

(17.97) 

3800 

(16.90) 

.15 

Three- in-series 

16.2 

3470 

(15.44) 

3360 

(14.95) 

.06 

The  Kolmogorov-Smlmov  goodness  of  fit  test  was  used  to  see  how  well 
the  Weibull  distributions  described  the  data.  In  all  cases  the  fitted 
distributions  were  acceptable  at  a  significance  level  of  5%. 

Fig.  4  presents  the  Weibull  CDF  for  the  fatigue  life.  The  solid 
curve  represents  the  CDF  for  the  basic  element  and  the  dashed  curve  the  three' 
in-series.  The  ten  data  points  for  the  basic  element  are  plotted  as  x's 
and  the  nine  data  points  for  the  three- in-series  are  given  as  open  circles. 

It  can  be  seen  that  the  scatter  for  both  the  basic  element  and  three-ln- 
series  specimen  is  large.  The  Weibull  distributions  fit  the  data  points 
fairly  well. 

Fig.  5  shows  the  CDFs  and  data  points  for  the  static  strength.  Once 
again  the  solid  curve  and  x's  are  for  the  basic  element;  the  dashed  curve 
and  open  circles  give  the  three-in-series  information.  The  static  strength 

is  seen  to  have  much  less  scatter  than  the  fatigue  life. 
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Figure  4  Distribution  Curves  for  Compression  Fatigue 
Data  of  Graphlte/Epoxy  Max.  Compressive  Load 
“  74%  Basic  Element  Mean  Strength 
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Figure  5  Distribution  Curves  for  Static  Conqpresslon  Strength 
Data  of  Graphlte/Epoxy. 

Q,  Basic  Element  Scale  Parameter  »  4040  Ibf  (17.97  kM) 
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Looking  at  the  estimated  population  means  we  find  that  the  static 
strength  of  the  three- in-series  specimen  is  89%  of  that  of  the  basic  ele¬ 
ment.  However,  the  fatigue  life  of  the  three-in-series  is  only  31%  of 
the  basic  element.  These  percentages  are  seen  to  be  in  general  agreement 
with  those  discussed  in  the  previous  section  that  were  calculated  from 
the  sample  means.  j. 


JL 
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VI.  Analysis  and  Conclusions 

From  the  theoretical  analysis,  It  can  be  seen  that  once  the  distribu¬ 
tion  of  the  basic  element,  Eq.  (2)  Is  known,  the  distribution  of  the  In- 
serles  specimen  can  be  calculated  from  Eq.  (3) .  This  Is  done  for  both  the 
static  strength  and  the  fatigue  cases,  and  the  results  are  presented  In 
Table  V.  They  are  also  shown  In  graphical  form  In  Figs .  6  and  7 •  The 
calculated  ones  are  labeled  "predicted",  while  the  fitted  experimental  ones 
are  labeled  "experimental". 


Table  V 

Comparison  of  Measured  and  Predicted  Properties 
of  Three  In-Serles  Model 


a.  Fatigue  Case  (Large  scatter,  small  shape  parameter). 

Three  In  series. 


Basic  Element  — 

Predicted 

Measured 

Shape  Parameter 

1.36 

1.36 

1.09 

Mean  Life 
cycles 

80400 

35900 

24700 

b.  Static  Case 

(Small  scatter,  large  shape  parameter) 

Three  In 

Basic  Element 

Predicted 

series 

Measured 

Shape  Parameter 

7.6 

7.6 

16.2 

Mean  Strength 
Ibf  (kN) 

3800 

(16.90)  • 

3290 

(14.63) 

3360 

(14.95) 
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Figure  6  Predicted  and  Measured  Distribution  Curves 
for  Compression  Fatigue  of  Graphite/Epoxy. 


Figure  7  Predicted  and  Measured  Distribution  Curves  for 
Static  Compression  Strength  of  Graphite/Epoxy 
6,  Basic  Element  Scale  Parameter  =  4040  Ibf  (17.97  kN) 
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For  the  life  of  the  three-in-series  specimen  we  see  from  Fig.  6  that  the 
shapes  of  the  experimental  curve  and  predicted  one  are  close;  however  the 
predicted  mean  life  is  45%  higher  than  the  experimental.  One  possible 
explanation  for  this  is  that  the  stress  distribution  is  non-uniform  along 
the  length  of  the  specimen  and  the  bottom  hole  is  subjected  to  higher  load, 
and  therefore  has  shorter  life  than  the  other  two.  This  is  born  out  by  the 
fact  that  in  the  three-in-series  fatigue  case  over  50%  of  the  specimens  failed 
at  the  bottom  hole. 

In  the  case  of  the  static  strength  the  shapes  of  the  experimental  and 
predicted  are  different,  but  the  mean  strengths  agree  to  within  3%. 

Based  on  the  results  and  analysis  presented  in  this  report  the  follow¬ 
ing  conclusions  can  be  drawn, 

A.  The  statistical  scatter  of  strength  and  life  of  a  basic  element  is 
one  of  the  main  reasons  for  the  decrease  in  strength  and  life  for  larger 
structures.  In  other  words,  the  scale  effect  in  fatigue  can  be  explained  by 
statistical  considerations. 

B.  In  general,  a  large  structure  contains  many  small  elements  which  are 
statistically  arranged  in  a  combination  of  in-series  and  in-parallel  modes. 

For  certain  special  structures,  a  simple  in-series  arrangement  can  be  assximed. 
In  the  present  case,  the  three-holed  specimen  is  truly  an  In-series  model, 
and  its  life  and  strength  calculated  from  theory  are  in  good  agreement  with 
measured  data. 

C.  The  larger  the  scatter  of  the  strength  or  life  among  specimens,  or  the 
smaller  the  Weibull  shape  parameter,  the  larger  their  decrease  for  large 
structures  due  to  the  scale  effect.  Since  in  general  there  is  more  scatter 
in  fatigue  life  than  in  static  strength,  the  scale  effect  is  more  pronotmcrd 
in  fatigue.  It  is  also  observed  that  there  is  more  scatter  in  fatigue  life 
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for  graphite  composite  materials  than  for  metals.  Consequently,  fatigue 
data  obtained  from  small  composite  coupons  must  be  used  with  caution 
when  applied  to  larger  structures. 


( 
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Nomenclature 


CDF 

C.O.V. 

f(x) 

F 

F 

F(x) 

n 


*0 

X 


X 


Xn 

x(n) 


a 

6 

0 

y 


Abbreviation  standing  for  cumulative  distribution  function 

Coefficient  of  Variation 

Probability  density  function 

Force  applied  to  a  bundle  of  fibers 

Breaking  force  of  a  bundle  of  fibers 

Cumulative  distribution  function 

Nximber  of  fibers  or  elements  (as  a  subscript  -  the  number  of 
elements  in  series;  as  a  super  script  in  parenthesis  -  the 
number  of  elements  or  fibers  in  parallel) . 

Also  total  number  of  data  points  in  a  sample. 

Number  of  surviving  fibers  or  elements 

Cumulative  distribution  function 

Median  rank 

Sample  standard  deviation 
Value  of  strength  or  life 
Position  parameter 
Sample  mean 

Value  of  breaking  force  of  a  fiber 
Strength  or  life  of  a  basic  element 
Strength  or  life  of  an  n  in-series  element 

The  failure  load  of  the  bundle/n. 

Weibull  shape  parameter 
Weibull  scale  parameter 
Population  standard  deviation 
Population  mean 
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Appendix  A  The  In-Parallel  Model  -  Bundle  Theory 

The  in-parallel  model  of  scale  effect  considers  a  number  of  basic  ele¬ 
ments  (fibers)  arranged  in  parallel  so  that  if  one  or  more  break,  the  sur¬ 
vivors  will  carry  a  redistributed  load.  The  theory  for  this  model  is  based 
on  the  bundle  theory  of  Daniels  [15].  The  model  is  applicable  to  both  the 
static  and  fatigue  modes  of  loading.  However  the  existing  set  of  equations 
of  the  bundle  theory  can  be  applied  only  to  the  static  case  and  is  reported 
in  this  appendix.  In  order  to  apply  the  bundle  theory  to  the  fatigue  case, 
additional  development  is  needed. 

The  Exact  Formula  The  cumulative  distribution  function  (CDF)  of  the  fiber 
under  a  static  load  is 

F^(x)  =  P(X  £  x)  (A-1) 


where  X  is  the  breaking  load  of  the  fiber.  Grouping  n  fibers  together  in 
parallel  will  form  a  bundle.  If  this  bundle  is  subjected  to  a  static  loading 
and  if  after  any  fiber  failure  the  survivors  share  the  load  evenly  then  the 
distribution  of  the  bundle  is 


x(n) 


P(X 


(n) 


n 


ix)  =  I 

m=l 


r 


I  (-l)"'“n!iP(X 


-ir. 


—  r 


1-1 


-  ri+r2/J 


[p(X  <_  x)J 


r*  f  •  •  < 
2' 


r  ! 
m 


(A-2) 


where  X^*^^  =  the  failure  load  of  the  bundle/n. 
greater  than  or  equal  to  one  such  that 

m 


The  t^'s  are  integers 


where  0  <  m  <  n.  The  inner  sum  in  Eq.  (A-2)  indicates  the  sum  over  all 


(A-3) 
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combinations  of  subject  to  the  condition  given  In  Eq.  (A-3) . 

For  two  fibers  arranged  In  parallel  Eq.  (A-1)  reduces  to 

j<  x)  =  2P(X  £  2x)P(X  £  x)  -  [P(X  £x)]^  (A-.4) 


As  the  number  of  elements  or  fibers  Increases  Daniels’  exact  formula  quickly 
becomes  unmanageable.  We  have  expanded  Eq.  (A-2)  for  f^alues  of  n  up  to 
n=7.  For  the  case  of  n=7  the  expression  contains  64  terms.  The  expanded 
formulas  and  the  associated  computer  program  listing  are  given  at  the  end 
of  this  appendix.  However,  as  the  number  of  elements  becomes  large, 

PCX'"  ^  x)  approaches  a  normal  distribution  and  its  mean  can  be  expressed 
in  a  simple  formula.  This  will  be  derived  below. 

The  Large  Bundle  Theory  The  large  bundle  theory,  which  was  presented  by 
Daniels,  can  be  derived  by  simple  considerations.  The  result  is  limited, 
however,  to  the  mean  strength  of  the  bundle.  The  standard  deviation  of  the 
normally  distributed  strength  was  given  by  Daniels  after  a  lengthy  statistical 
derivation.  In  the  following,  we  shall  derive  the  mean  strength  expression 
and  present  Daniels’  standard  deviation  formula  without  derivation. 

Consider  a  bundle  of  fibers  subjected  to  a  force  F.  Each  fiber  will 
experience  a  force  x  where 


X 


(A-5) 


and  n’  is  the  total  number  of  surviving  fibers  currently  in  the  bundle. 

Let  the  probability  that  a  single  fiber  breaks  under  the  force  x  be 
P(X1  x). 


P(X  <  x) 


'X 
-  0 


f(C)d? 


(A-6) 


where  f(t)  is  the  probability  density  function. 
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The  reliability  is  given  by 


R(x)  =  1  -  P(X 


*  "vr 


f(?)d5 


The  total  niomber  of  fibers  surviving  will  be  n' , 
n'  =  n  R(x)  «  n[l  -  P(X  £x)] 


(A-7) 


(A-8) 


where  n  is  the  total  number  of  fibers  originally  in  the  bundle. 

The  relationship  between  the  force  applied  to  the  bundle,  the  force 
experienced  by  each  surviving  fiber,  and  the  number  of  original  fibers  is 
obtained  by  substituting  Eq.  (A-8)  into  Eq.  (A-5)  and  rearranging,  or, 

F  =  n  x[l  -  P(X  x)]  (A-9) 

The  breaking  load  can  be  obtained  by  maximizing  the  force  F  in 
Eq.  (A-9) ;  this  yields 

dF 


dx 


=  0 


or 


dx 


{nx[l  -  P(X  £  x)]}  =  0 


(A-10) 


Performing  the  differentiation  yields  the  equation 

[1  -  P(X  <_  x)]  =  X  f(x)  (A-11) 

The  value  of  x  that  solves  Eq.  (A-11)  will  be  designated  as  x^.  The 
'corresponding  value  of  the  breaking  load  is  obtained  by  substituting  x^ 
into  Eq.  (A-9) .  Thus 

F^  =  n  x^[l  -  P(X  £  x^)]  (A-12) 

Daniels  shows  that  large  bundles  have  failure  strength  that  are  of 
a  normal  (or  standard)  distribution.  The  mean  strength  is  then  given 

by  F^/n,  and  the  standard  deviation  is 
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,(n) 


.^■(X  £  x^)  [1  -  P(X  < 


TJT 


(A-13) 

When  a  Welbull  distribution  is  assumed  for  the  single  fiber,  the  value 
of  can  readily  be  obtained  from  Eq.  (A-11) .  The  Weibull  density  and 
distribution  functions  are 


f(x)  = 


a(x-x-)' 


,a-l 


exp 


and 


P(X  ^  x)  =  1  -  exp 


;C-?) 

■  C-?)' 


(A-14) 


(A-15) 


We  will  further  assume  x^  =  0. 


Substituting  the  Weibull  distribution  into  Eq.  (A-11)  and  solving 
gives  x^ 

J _ 


X  = 

T 


(ct) 


1/a 


(A-16) 


The  breaking  force  from  Eq.  (A-12)  will  therefore  be 


F  -  n  ^ 

-  n  - rT~  exp 


(a) 


1/a 


(0 


(A-17) 


Then  for  large  bundles,  the  mean  and  standard  deviation  will  be 

,.(n) _ e 


(a) 


(A-18) 


and 


,(n)  _  3 


1/a 


[1  -  exp(-  i)]  exp(-  i) 


(A-19) 
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k 


r 


The  coefficient  of  variation,  C.O.V. 


is  just  the  ratio  of  0  to  y.  Therefore 


for  a  large  bundle 


C.O.V. 


(n)  _ 


[1  -  exp(-  — )] 


n  exp(- 


1/2 


(A-20) 


As  can  be  seen  the  mean  strength  is  independent  of  the  number  of  fibers 
(or  olomento)  arranged  in  parallel  whereas  and  C.O.V. are  Inversely 

proportional  to  the  square  root  of  n. 

Figures  A-1-3  show  the  CDF  of  the  bundle  where  the  basic  elements  are 

Weibull  with  indicated  shape  parameter.  The  distributions  for  n  =  2,3,4, 
5,6,7,50,100,  and  1000  elements  in  parallel  plus  the  basic  element  are 
shown.  Values  of  shape  parameter  of  the  basic  element  used  are  2,14,  and 
24.  The  plots  have  been  nondimensionalized  by  the  scale  parameter  of  the 
basic  element.  Also  tabulated  on  these  figures  are  the  values  of  y/6,  o/P. 
and  C.O.V.  for  each  arrangement. 

Interpolation  for  n=7  to  50  For  small  values  of  n,  up  to  n=7,  we  have 
expanded  and  calculated  the  exact  bundle  formula.  For  large  values  of  n, 
the  large  bundle  equations  are  applicable.  There  is  a  need  to  have  formulas 
for  the  intermediate  range  of  n,  say,  between  n=7  and  50.  We  shall  attempt 
to  do  this  by  interpolation.  Figure  A-4  is  a  plot  of  the  mean  strength  for 
both  large  bundles  and  small  bundles.  The  values  for  n=7  to  50  can  be 
interpolated  as  shown.  Similarly,  Fig.  A-5  shows  the  corresponding  plots 
for  the  coefficient  of  variation. 

PiannPd  Experimental  Work  In  our  experimental  program,  we  intend  to  test 
Daniels'  formula  for  three  elements  in  parallel  subjected  to  a  static  load. 
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Distribution  Curves  of  n  Basic  Elements  in  Parallel 


Cumulative 

Distribution 


Distribution  Curves  of  n  Basic  Elements  in  Parallel 


Strength 
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Figure  A-4  Mean  Strength  for  n  Elements  in  Parallel 


Coeficient  of  Variation  for  n  Elements  in  Parallel 
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The  in-parallel  model  subjected  to  a  fatigue  loading  is  considerably 
more  involved  than  the  static  case.  This  is  primarily  due  to  the  fact 
that  after  some  fiber  failures,  the  loading  in  the  remaining  fibers  in¬ 
creases,  and  the  corresponding  life  distribution  also  changes.  Further 
analytical  work  must  be  done  before  an  expression  to  predict  fatigue  life 
can  be  presented.  The  experimental  portion  of  this  phase  presents  some 
special  problems.  One  problem  is  the  transient  loading  following  a  fiber 
failure  may  cause  additional  failures.  Another  problem  is  that  the  load 
in  the  surviving  fibers,  for  n=3  or  4,  may  exceed  the  static  strength.  One 
testing  procedure  which  may  prove  successful  is  to  operate  in  the  constant 
stroke  mode.  The  minimum  applied  load  will  be  10%  of  the  maximum.  The  three 
elements  will  be  placed  side  by  side  in  the  gripping  anti-buckling  apparatus. 
When  one  specimen  fails  the  test  will  be  stopped  and  the  remaining  specimens 
rearranged  so  as  to  insure  uniform  loading.  The  test  will  then  be  restarted 
and  run  until  the  next  specimen  fails.  After  positioning  the  last  element 
the  test  will  proceed  to  total  failure.  The  total  life  of  the  three-ln-parallel 
arrangement  will  be  recorded. 

Expanded  Formulas  for  Small  Values  of  n 

Eq.  (A-2)  has  been  expanded  for  values  of  n  up  to  7.  On  the  following 
pages  these  expansions  are  given.  Following  these  is  the  computer  program 
listing  used  to  evaluate  ^  x)  ,  ^  and  ^  • 
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E  X  pclrsion  O 

for  n  =•  v3 


^  C>CPbCP^^Pt^  —ZCP/^CPif  —  30/>V/='3> 

Py  “  P^X  ^  v)  "  /  —  ^ 

P2=  por^  ^)  =  /  -  ^yo£-i^-^Sj 

P3=  p<2f>  3x)=  /  -  ^/- £' 
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Exp««S»n  of 
•for  A  =  ^ 

p/  =r  /=^x£.x>  =  ^  “  J 

P2,= p^xi6.^  ^  - 1  - 

pz= pac±^^'^-i 

r>-^  -  p<ix  i  — I  —  <270  £'-<r'^^^)  J 
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df  E<^.  (A^z) 

h=  v5“ 


K)  =■  -  /oCP4f>^^^d^ 

““  ^oCp^CPi^  -f-  ZoCPS)Cp^)CPd^ 2a 

Zo  CPt^CPt:xlPh  '^3o  Cp^^CP^CPd'p^ocp^'9‘CP^P^y 

^  'SoCPS'^CPtf^CPif^  -<B06PS^CP^:)<P2^CPy'^^ 

p 

-&o  CPs^yi)CP:^Cpy')  ^Psr^O>^^  <Pz^cpJ^ 

— 6o  Cp^^Cp^Cp^  CPi^-i-  J2jo  dPS>CP4y>^:^z:^CP'^ 

p>/  -  -?=C5C  ^  =^  y  -  ty^€-<^^T  1 

PZ--  POC£  %?■')=  I  - 
pz=  p^ic  £  %)  =/ 
p^ = pac£  %) 

PS-  =  P<^£~S->^'^  -/  - 
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E/payiswi^ 

for 


O’®  ^ 


^^oCf^"^CP^(>*)? -  ^CP^J'CPt^CPi^ 

-t  •l-l2o0=*^)O*S^t}CP*^  tTJDCPC^P^iCP^^*^ 

-f'/2o^:^^>73y/>^^)  ^SoCPS^'^Cf^^^CP/^ 

'+‘iSoCP^O»t^0^i^Pb  -r  /SoCP^PSrXP^CP'f' 
+i^oCPC)O^CPt}0^h  'h  tBoO^XP^CP^C^^^ 

-  ■s<i  O  O^^^CP4X^i^CP  25 
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Pj  =  PCX  ^^'>  =  /  ~  ey,£’-^'Tj 

Pz  =  PCX :£  %>  =  >  -e.^r-(P^pT'j 

P3-  PCXC 

P'^^PCX  ~  —i~  y 

ps=j>cx^  ‘^>=/  -  ^L-c-tjpy'j 

p<o^  PCX^  C>)c^  —  /  — 
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ExpiS»>s iai' 

for  A  =-7 

^p>')’  -9^i='?'^c^''f  -oCPzo^<:p''> 

-  Z.I CPS'f/lP'f  ~^CP0'>^CPd^ 

V-  ^zO‘'/^CP^'><>>f->-  ^zCP'y')0>^<Pf>  +  ^zCPtfCPioCPy'o 
-^>os0>i'^-0‘s-yp>;f  -h  ic^o>^y:psy<pf^->-/'=ccp<it<:^ii‘'cp''> 
~PJ^O>P')CPZ'^‘'Cph^  ■h/osCf^'fo>7'i^<^Ph  -t-Jos<;p4o‘‘ CPtoCP-/: 
+  2J  oaxi^Cp^o/f  ■pz^c>(>&fcpsfzp'i'*z'o<:psf<p-i?o:p/o^ 

•h  -p/^o<CP'>yp*f<lP'^ 

-z,oo>7yp6),0‘£^<^‘^'^  cpo-^t>&ypz^‘><:pd 

-no  cp^yp^'^o’tyjp^  ~  ^<5  cf^cptocpto  cppi 
-i4  zjc>  cp7yipc,ypp>‘^^p^^  —  ^  ^0  o>7y/>c.yp^  o^i)^ 

-XJZ.O  CP'^O‘S'0^  —^toO>C^CF^O:‘^hCP''^ 

-nnoCPifCP^^CJ'iiCPh  z^CP-h O4o^ Cpi^Lp^ 

—4/2-0  0=<hCP^'^CP^^P^  -*tzt>  i!.PSr:^Cp^&'nCPi'> 
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—^ZaCf^S':^CP^'^O^t^CPi^  --4'Z'O 

k 

-f-^o  CP*bCF>(^O^s^0^4)Cp^^  v-  B^/o 

■r^^o  Cf*^yO=>f^ 

-2. 

•f*  V*  j  ?<so 

-hneo  j  z<^o 

'hi  -zGo ■^jZ(ao  o=>(^Y</^s>0*^^2:>o^d 

-  2S2-oCp^?^(^0=^^Cp'^0*t^O^^'> 

-zB2  0C^r^^o>^y/=^^yp^:)O^i^  -zszc 

'^ 5^04^0  O^j) 
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7^/ = 

-  ;>CX£  z  j  -  J 

¥)  =  /-  '^fJ 

p>-^-  p  cx  ^  /  -  ^.y4^-<r  ^/sy^j 

PS  -  PCX  y  I  - 

^  P^.X  i  -7  -  1 

/=7  =  P<:X  ^  '^y'y  ~/  —  J 
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Computer  Listing  for  ^ 


and  0 


(n) 


OPfinPA  /MF7336 


start  of  segment 

-> - —  OTMpTi s  TfTR  Ti  raoo T  T T 3(nn70  r3 0 rtT  w  c  30 0 ~3  fiTv 1 3 ffo-TTvv  c-jotn >' - 

IVANC^CO  >  ^ 

- NP-»-m - : - 

■>  BaifA. 

- - - - - - - 


Fortran  for  expansion  of  Eq,  (A-2)  for  n  =  2,  3,  4, 
5f  6,  7  Is  given  on  the  following  pages  and  Is  to  be 
Inserted  In  this  location  of  program. 


> 


;2- 


>  ■ 

- 88 - CONTTNUE - - 

50  I»l 

>  5t - coNritinf - - - - - 

IFtT.LT.3)Gn  TO  1000 

— - rr(TirGTVTnTGcrTcr'7ot)o - 

y  IFCt.6T.NP-5)Gn  TO  1500 

- XT“0<T-7> - - - - - - - - - - - - ^ - 

>  ■: - X-3»otrj - - - - - -■■■■■■- - - 

X4a0fl+1) 

■  ^ ^ - X5-o<-f-»73 — ^ ^ ^ ^ - : — - — ^ — — - — - ^ - 

>  X6®o(T*3) 

- iriaTrr-7T - - - 

Y?»v<i-n 

y - — - YSaY-rtT - - — - - ^ - — - — - 

YAaYfl+l) 

- Y5»ytt>1D - 

3  Y6»V(I*3) 

- 69—  XsXi - - - - - ; - 

YA*V3 

^  - - - IFtT.FO.NP-3)Gn  T(r-2000 -  -  — - - 

100  TFKm18((X-X?)*(X-X4)*(X-X5)*(X-XA)*Y1 1/((Y1»X2)*C vi-X3)*(Xl-X4)*( 

- tX  1»  Xx;  )-«  (  X  l'Xiy)~Y - 

P  110  TEPM9s( (X-Xn*(X-X4)*(X-X5)*(X-V6)*Y7)/( (X9-X1 ) * ( X7-X3 ) * < X7-X4 ) ‘(X 

17-y7v*CX?*X53  J -  - - 

120  0FN-»r(X3-Xn*CX3*X2)*(X3-X4)*CX3-X5)*tX3-y6) 

t.-  - - 121  -  TFKM3A  =  (y-X7)»(X-X«l  )*tX-X5)*(  X-Xft)*Y3/0FN3 - - 

'  122  TFKM3Raf X-Xt 5*(X-X4)*(X-Xb)»{X-yA)*Y3/0rN3 

- 123 - TEl^H3C  a(  X-XT)TtX  X-X2  )  «T  X-*53*t  X-X  6  )*Y  3-/eFrMy - 


P  124  TFKMSDaCX'XI )*(X-x2)*(X-X«)*(X-yfi)*Y3/orN3 

- 125  TEf'M3F  =  (  X-X1  )  ♦  (  X-X2  )  *(  X-X4  )  •  t  X-X5  )*  Y 3/OF N  V 

127  TFKM3sTFR-l34  +  TFRIi3S+rERM3C  +  TERM30  +  TFRM3F 
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130 
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- TF  K  H2  =  C  X  -  X  3  >  *  r  X  “■  X  4  )  » t  X  •XS)*tX“Y  fi  1  *  Y?  /0FN2- 
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Fortran  for  In-Parallel  Model,  n  =  2 
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Fortran  for  In-Parallel  Model,  n  =  3 
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Fortran  for  In-Parallel  Model,  n  =  5 
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Appendix  B 


Detail  Drawings  for  Gripping  Anti-Buckling  Guide  Device 


55 


'igure  B-2  Rear -Grip 

Dimensions  in  Inch,  1  in  =  25,4  mm) 
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Hat  er lal :  Aluminum 


Figure  B-3  Front  Anti-Buckling  Guide 
(All  Dimensions  in  Inch,l  in  =  25.4  mm) 


Material:  Aluminum 


Figure  B~4  Rear  Anti-«Buckling  Guide 
(All  Dimensions  in  Inch,  1  =  in  =  25.4  mm) 
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Appendix  C.  Additional  Data 

Some  additional  fatigue  tests  were  performed  with  elements  having 
different  geometries.  All  tests  were  performed  with  the  same  loading  con¬ 
ditions  as  described  in  the  main  body  of  this  report.  The  first  set  of  tests 
was  performed  with  elements  containing  four  1/4  inch  (6.4  mm)  holes  as  shown 
in  Fig.  C-1.  Table  C-I  gives  the  life  and  the  failure  location  of  these 
elements.  The  mean  life  is  5430  cycles.  This  is  much  less  than  the  in¬ 
series  model  prediction  life.  Due  to  the  fact  that  all  elements  failed 
through  the  bottom  hole  it  was  concluded  that  its  location  was  too  near  to 
the  end  tab  and  thus  experienced  a  greater  load  than  the  other  three. 

Several  months  later  a  series  of  fatigue  tests  were  performed  on  ele¬ 
ments  containing  two  1/4  inch  (6.4  mm)  holes  as  shown  in  Fig.  C-2.  The 
initial  tests  in  this  series  were  performed  with  configuration  1  type 
elements.  These  all  broke  through  the  bottom  hole.  The  geometry  was 
altered  by  moving  the  two  holes  further  apart  as  shown  as  configuration  2. 
These  also  all  broke  through  the  bottom  hole  and  had  shorter  life  than 
configuration  1  elements.  The  lives  for  these  two  configurations  are  shown 
in  Table  C-II, 

In  order  to  see  if  the  position  of  the  hole  in  an  element  affected  the 
life,  a  series  of  three  tests  were  performed  with  an  element  with  a  1/4 
inch  (6.4  mm)  hole  located  near  the  bottom  end  tab  as  seen  in  Fig.  C-3. 
l3-ble  C— III  records  the  life  for  these  tests.  As  can  be  seen  the  lives  are 
much  shorter  than  the  lives  recorded  for  the  basic  element  (Table  la). 


Table 
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(All  Dimensions  in  Inch,  1  in  ==  25.4  mm) 
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(All  Dimensions  in  Inch,  1  in  =  25.4  mm) 


Table  C-III 
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Figure  C-3  Single  Hole  Element 
(All  Dimensions  in  Inch,  1  in  =  25.4  mm) 


WADC.7«259^0 


This  tends  to  indicate  that  the  closer  the  hole  to  the  bottom  end  tab,  the 
greater  will  be  the  force  seen  by  the  hole  and  hence  the  shorter  the 
life. 

The  results  obtained  in  Tables  C-I-III  tend  to  indicate  that  the 
stress  distribution  is  not  uniform  throughout  the  total  length  of  the 
specimen I  Since  the  dynamic  loading  is  applied  by  the  accuator  to  the 
bottom  of  the  specimen,  the  stress  concentration  at  the  bottom  hole  is 
greater  than  at  the  other  holes. 

This  evidence  is  of  course  not  conclusive.  Other  plausible  explana¬ 
tions  can  be  forwarded.  First  the  original  basic  and  three-ln-serles 
elements  were  taken  from  one  end  of  the  plate  whereas  the  elements  discussed 
in  this  paragraph  were  obtained  from  the  other  end  of  the  plate.  It  is  ^ 

conceivable  that  the  two  ends  of  the  plate  had  different  properties  and 
thus  the  elements  taken  from  each  constituted  distinct  populations.  A 
randomized  scheme  of  picking  test  elements  of  course  would  have  eliminated 
this  possibility.  Another  explanation  was  that  the  tests  listed  in 
Tables  I,  II  and  Tables  C-I-III  were  performed  before  and  after  respectively 
a  major  breakdown  and  repair  of  the  Instron  testing  machine.  Also,  the 
gripping  anti-buckling  guide  may  have  been  inserted  in  the  Instron  improperly.  * 

In  short,  due  to  the  uncertainties  associated  with  Tables  C-I-III  no  further 
analysis  was  attempted  on  this  data.  ** 


64 


N  ADC -78259 -60 


DISTRIBUTION  LIST 
Government  Activities 


Ko.  of 
Conies 


pf  a!  Parmley)  .  .  •  .  - - -  ^  ^  ; 

(Attn:  FBC/Mr.  C.  Wallace)  . . . 

(Attn:  FBC/Mr.  E.  E.  Zink)  . 

AFML,  WPAFB,  OH  A5433 

(Attn:  LAM  (Technical  Library)).  ...•••• 

(Attn:  LT-l/Mr.  W.  R.  Johnston) . . . .  * 

(Attn:  LTF/Mr.  T.  Cordell)  . . !  .  .  * 

(Attn:  FBSC/Mr.  L.  Kelly). 

(Attn:  MAC/Mr.  G.  P.  Peterson)  . 

(Attn:  MXA/Mr.  F.  J.  Fechek)  .  .  . . 'll.. 

(Attn:  MBC/Mr.  T.  G.  Rcinhard,  Jr.) . 

AFOSR,  Washington,  D.C.  20333  ^  . 

(Attn:  Mr.  J.  Pooerantz) . '  1  1  1  1  1  1  1  I  •  •  • 

FAA*  Airfraies  Branch’.  ES-12ol  Washington,  D.C.  20^3 

(Attn:  Mr.  J.  Dougherty)  •  •  ;  •  ..•••••••••••• 

NASA  (ADM),  Washington,  D*C.  20546  ^  ^  ^  ^  ^  ^ 

NASAt*^GeSrge*c!^S2hSll  Space  klght  Center,  Huntsville, 

AL  35812  .  ^  - 

(Attn:  S&E-ASTO-ES/Mr.  E.  E.  Engler)  .....  . 

(Attn:  S&E-ASTN-M/Mr.  R.  Schwinghamer) . .  .  .  .  • 

(Attn:  S6£-ASTM-MNM/Dr.  J.  M.  Stuckey)  •  •  •  •  ••••••• 

NASA,  Lewis  Research  Center,  Cleveland,  OH  44U  _ 

(Attn:  Technical  Library,  and  M.  Hershberg).  ....•• 
NAVPGSCHL,  Monterey,  CA  95940  .  .  ,  x  .  .  .  . 

(Attn:  Prof.  R*  Ball,  Prof.  M.  H.  Batik) . .  • 

NAVSEASYSCOM,  Washington,  D.C.  20362 

(Attn:  Code  035,  Mr.  C.  Pohler) . . 

NAVSEC,  Hyattsvllle,  MD  20782 

(Attn:  Code  6101E03t  W.  Graner)# 

NAVSHIPRANDCEN,  Bethesda,  MD  20034 

(Attn:  Code  173.2,  Mr.  W.  P.  Cauch) . 

NAVSHIPRANDCEN,  Annapolis,  MD  21402  ^  ^  ^ 

(Attn;  Code  2870,  Mr.  H.  Edelstein).  .  . . .  •  • 

NOL.  White  Oak,  MD  20910  _  ^ 

(Attn:  Mr.  F.  R.  Barnet) . . . 

NKL,  Washington,  D.C,  20375  .  .  •  .  . 

(Attn:  Dr.  I.  Wolock). 

CNR,  Washington,  D.C,  20362 

(Attn:  Dr.  N.  Perrone)  . . . . 


NADC-78259.60 


Government  Activities  (Cont.) 

FIASTEC,  Plcatinny  Arsenal,  Dover,  HI.  07801 

(Attn:  Librarian,  Bldg.  176,  SARPA-FR-M-D  and  Mr.  H.  Pebly).  .  2 

Scientific  &  Technical  Information  Facility,  College  Park,  MD 


(Attn:  KASA  Representative).  •*••••••••••••••••  1 

OSAAVMATLAB,  Fort  Bust is,  VA  23603 

(Attn:  Mr.  R,  Beresford)  1 

USAMATRESAG,  Watertown,  MA  . 

(Attn:  Dr.  E.  Lenoe)  1 

DSARESOFC,  Durham,  NC  27701  . . 1 


Mon-Govemnent  Agencies 


Avco'Aero  Structures  Division,  Nashville,  TN  37202 

(Attn:  Mr.  W.  Ottenville).  ..................  1 

Battelle  Columbus  Laboratories,  Metals  and  Ceramics  Information 

Center,  505  King  Avenue,  OH  43201 . 1 

Bell  Aerospace  Company,  Buffalo,  NY  14240 

(Attn:  Zone  1-85,  Mr.  F.  M.  Anthony)  .............  1 

Bell  Helicopter  Company,  Fort  Worth,  TX  76100 

(Attn:  Mr.  Charles  Harvey)  .................. 1 

Bendix  Products  Aerospace  Division,  South  Bend,  IN  46619 

(Attn:  Mr.  R.  V.  Cervelli)  . . .  1 

Boeing  Aerospace  Company,  P.O.  Box  3999,  Seattle,  WA  98124 

(Attn:  Code  206,  Mr.  R.  E,  Horton)  ..............  1 

Boeing  Company,  Renton,  Washington  98055 

(Attn:  Dr,  R.  June).  ..  .  1 

Boeing  Company,  Vertol  Division,  Phila. ,  PA  19142 

(Attn:  Mr,  R.  L.  Pinckney,  Mr,  D.  Hoffs tedt)  .........  2 

Boeing  Coiqiany,  Wichita,  KS  67210 

(Attn:  Mr.  V.  Reneau/M5  16-39)  1 

Cabot  Corporation,  Billerica  Research  Center,  Billerica,  MA 

01821  . . . . .  1 

Drexel  University,  Phila. ,  PA  19104 

(Attn:  Dr.  P,  C,  Chou)  .....  .  1 

E.I.  DuPont  Company,  Wilmington,  DE  19898 

(Attn:  Dr.  Carl  Zweben)  Bldg.  262/Room  316  1 

Fairchild  Industries,  Hagerstown,  MD  21740  “  t.—  ——— 

(Attn:  Mr,  D,  Ruck).  1 

Georgia  Institute  of  Technology,  Atlanta,  GA 

(Attn:  Prof,  W,  H.  Horton)  . .  1 

General  Dynamics /Convair,  San  Diego,  CA  92138 

(Attn:  Mr.  D,  R.  Dunbar,  W,  G.  Scheck)  . . .  2 

General  Dynamics,  Fort  Worth,  TX  76101 

(Attn:  Mr.  P.  D.  Shockey,  Dept.  23,  Mall  Zone  P-46) . 1 

General  Electric  Company,  Phila.,  PA  19101 

(Attn:  Mr.  L.  McCreight)  . . .  1 

Great  Lakes  Carbon  Corp.,  N.Y.,  NY  10017 

(Attn:  Mr.  W,  R,  Benn,  Ifer. ,  Markey  Development)  1 

Grumman  Aerospace  Corporation,  Bethpage,  L.I,,  NY  11714 

(Attn:  Mr,  R.  Hadcock,  Mr.  S.  Dastin) . . . 2 
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Kon'Govemnient  Agencies  (Cent.) 

Hercules  Powder  Company,  Inc.,  Cumberland,  MD  21501 
(Attn:  Mr.  D.  Hug)  . . . 

H.  I.  Thompson  Fiber  Glass  Company,  Gardena,  CA  90249 
(Attn:  Mr.  N.  Myers) 

ITT  Research  Institute,  Chicago,  IL  60616 

(Attn:  Mr.  K.  Hofar)  . . . . 

J.  P.  Stevens  &  Co.,  Inc.,  N.Y.,  NY  10036 

(Attn:  Mr.  H.  I.  Shulock).  . . 

Kaman  Aircraft  Corporation,  Bloomfield,  CT  06002 
(Attn:  -  Tech.  Library). 

Lehigh  University,  Bethlehem,  PA  18015 

(Attn:  Dr.  G.  C.  Sih) . . 

Lockheed-Callfomia  Company,  Burbank,  CA  91520 
(Attn:  Mr.  E.  K.  Walker,  R.  L.  Vaughn) 

Lockheed-Georgia  Company,  Marietta,  GA 

(Attn:  Advanced  Composites  Information  Center,  Dept.  72-14. 

Zone  42).  . . . 

LTV  Aerospace  Corporation,  Dallas ,  IX  75222 

Ur.  0.  E.  Dhonau/2— 53442 ,  C.  R.  Foreman).  ..... 
Martin  Coii5)any,  Baltimore,  MD  21203 

(Attn:  Mr.  J.  E.  Pawken) . . 

Materials  Sciences  Corp.,  Blue  Bell,  PA  19422  ....!!  IT  i 
McDonnell  Douglas  Corporation,  St.  Louis,  MO  63166 

(Attn:  Mr,  R.  C.  Goran,  0.  B.  McBee,  C.  Stenberg).  .  .  .  .  < 
McDonnell  Douglas  Corporation,  Long  Beach,  CA  90801 

(Attn:  H.  C.  Schjulderup,  G.  Lehman)  . 

Minnesota  Mining  and  Manufacturing  Company,  St.  Paul,  MN  55104 
(Attn:  Mr.  W.  Davis) 

Northrop  Aircraft  Corp. ,  Norair  Div. ,  Hawthorne,  CA  90250 

(Attn:  Mr.  R.  D.  Hayes,  J.  V.  Noyes,  R.  C.  Isemann).  .  .  .  . 
Rockwell  International,  Columbus,  OH  43216 

(Attn:  Mr.  0.  G.  Acker,  K.  Clayton).  ............ 

Rockwell  International,  Los  Angeles,  CA  90053 
(Attn:  Dr.  L.  Lackman) 

Rockwell  International,  TUlsa,  OK  74151 

(Attn:  Mr.  E.  Sanders,  Mr.  J.  H.  Powell) . .  .  .  .  . 

Owens  Coming  Fiberglass,  Granville,  OH  43023 

(Attn:  Mr.  D.  Mettes) . . 

Rohr  Corporation,  Riverside,  CA  92503 

(Attn:  Dr.  F.  Riel  and  Mr.  R.  Elkin)  ....  . 

Ryan  Aeronautical  Company,  San  Diego,  CA  92112 

(Attn:  Mr.  R.  Long).  .  . . . . 

Sikorsky  Aircraft,  Stratford,  CT  06497  •  .  .  .  .  .  . 

(Attn:  Mr.  J.  Ray) 

University  of  Oklahoma,  Norman,  OK  93069 

(Attn:  Dr.  G.  M.  Nordby)  . . . 

Union  Carbide  Corporation,  Clevelaivl,  OH  44101 
(Attn:  Dr.  H.  F.  Volk)  . . 
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